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Scaling functions for nonequilibrium fluctuations: 
A picture gallery 

Zoltan Rac^ 
Institute for Theoretical Physics, Eotvos University, 
Pdzmdny setdny 1/a, 1117 Budapest, Hungary 

The emergence of non-gaussian distributions for macroscopic quantities in nonequilibrium steady 
states is discussed with emphasis on the effective criticahty and on the ensuing universality of dis- 
tribution functions. The following problems are treated in more detail: nonequilibrium interface 
fluctuations (the problem of upper critical dimension of the Kardar-Parisi-Zhang equation), rough- 
ness of signals displaying Gaussian 1/f power spectra (the relationship to extreme- value statistics), 
effects of boundary conditions (randomness of the digits of vr). 
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I. INTRODUCTION 



Spatially averaged (global) quantities in homogeneous equilibrium systems display Gaussian fluctuations since the 
correlation length is, in general, finite and the central limit theorem applies. The situation is more complicated 
^ I at critical points where the long-range correlations yield diverging fluctuations which, in turn, generate nontrivial 
probability distribution functions (PDF-s). A remarkable simplification occurs, however, even in critical systems 
where the large-scale fiuctuations lead to universality For the PDF-s of macroscopic quantities, this means that 

^ the shape of the PDF-s depends only on a few general characteristics of the system (dimension, symmetries, range 
\ of interactions) 0, 0| • The emergence of universality from large-scale fluctuations appears to be so robust that one 

■ expects that similar mechanisms works in far from equilibrium steady states as well. Since no general theory exists for 
O nonequilibrium systems, studying the differences between equilibrium and nonequilibrium phase transitions can indeed 
^ I be instrumental in understanding some distinguishing but still robust properties of nonequilibrium systems 0, 0, 01 ■ 

In particular, it may help making inroads in the largely unknown territory of nonequilibrium PDF-s. 

At first sight, universality ideas should be of restricted use since they apply only to critical points. One should 
^ remember, however, that nonequilibrium systems displaying power law behavior in their various characteristics (cor- 
, relation in space or time, fluctuation power spectra, size-distributions, etc.) are abundant in nature. Examples range 
from interface fiuctuations Q and dissipation in turbulent systems Q to voltage fiuctuations in resistors 0, and 
to the number of earthquakes vs. their magnitude The underlying reason for the ubiquity of power-laws is not 

■ understood (the widely used expression self- organized criticahty 12] is a testimony for this fact) but the observed 
' effective criticahty suggests that a classification of nonequilibrium PDF-s can be developed using the logics of critical 
, phenomena. Namely, strong fiuctuations and power-law correlations imply universal scaling functions for the PDF-s 

^ ' and, consequently, the nonequilibrium PDF-s in a large number of phenomena can be determined by studying the 
^ , nonequilibrium universality classes. 

^ ' Compared with equilibrium systems, complications are expected to arise from the fact that the properties of 
nonequilibrium steady states are determined not only by the interactions but by the dynamics as well. Thus, dynamical 
symmetries (conservation laws, the effects of breaking of time-reversal symmetry, etc.) should also play an important 
Q ' role in the classification. Furthermore, in spite of being universal, the scaling functions do depend on the boundary 
Q : conditions Although this is an extra complication, it indicates that the scaling functions may be suitable for 

^ ' describing an important feature of nonequilibrium states, namely, that the bulk behavior depends on the boundary 
. ^ , conditions (note that fluxes are often generated by an appropriate preparation of the boundaries). 
^ ■ General (field-theoretic) studies of nonequilibrium universality classes 0, 0, do not address the question of 
' distribution functions due to technical difficulties. For practical purposes, on the other hand, it is important to have 
- - 1 the scaling functions associated with the PDF-s since, as we shall see below, they provide a possibility for "fit-free" 
comparisons with experiments. Once a gallery of such scaling functions has been built, it can be used to identify 
symmetries and underlying dynamical mechanisms in experimental systems; to discover analogies between seemingly 
different systems due to both belonging to the same universality class, and the applications are restricted only by 
imagination. 
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The project of building the picture gahery of scahng functions associated with nonequihbrium PDF-s has been going 
on fperhaps unknowingly) for several years. For simple systems, the scaling functions have been found analytically [Wi 
IrR ll^ llfll l2(l I21L l2^. In most of the cases, however, the PDF of the "microscopic" configurations is unknown and 
thus one has to resort to other means of calculations. One of them is a phenomenological approach which 
consists of introducing an effective Gaussian action with singular dispersion and fixing the dispersion to yield the 
observed scale- invariant fluctuations. The other is the brute force simulations of models which are believed to be in 
the same universality class as the system at hand [2^ |2^ . The resulting picture gallery is far from complete but 
contains a few interesting pieces which will be presented below. First, I will show how the scaling functions are defined 
and calculated using simple examples from surface growth problems (Secnj- The details of the derivation will be 
demonstrated on the example of a one dimensional {d = 1) surface dynamics that is equivalent to the problem of 
Gaussian 1// noise (Sec lIII CI . As we shall see, this calculation establishes a connection between the the 1// noise 
and one of the limiting distributions of extreme statistics. Applications of the scaling functions will be discussed in 
Sec lIII B1 with details presented in connection with the problem of the upper critical dimension of the Kardar-Parisi- 
Zhang (KPZ) equation. Finally, as a demonstration of the importance of the boundary conditions, we shall discuss 
the problem of the randomness of the digits of tt fScc lIVIl . 



II. SURFACE GROWTH AND SCALING FUNCTIONS 



Among the nonequihbrium systems displaying "effective" criticality, growing^urfaces provide a conceptually simple 
and versatile laboratory from both experimental and theoretical point of viewj^ . The criticality here means that 
growing surfaces are usually rough i.e. the mean-square fluctuations of the interface diverge with system size. More 




FIG. 1: Surface formed by deposition and evaporation. The vertical velocity of the surface = dth{x,t) is assumed to depend 
on the local properties, Vz(dxh, d^h, ...) and, at a macroscopic level of description, only a low-order gradient expansion of Vz 
is kept. The growth models differ from each other by the physical processes which are considered to be relevant and thus by 
keeping only the appropriate terms in the gradient expansion. 



precisely, let us denote the height of the surface above a d dimensional substrate of characteristic linear dimension L 
by h{f,t) (see Fignj. Then the mean-square fluctuations are defined as 

r 

where the spatially averaged height is given hy h = ^ph{f, t)/AL with Al being the area of the substrate. For rough 
surfaces one finds that the steady state average of the fluctuations {w2)l diverges as {'W2) l ^ L,^^, and the critical 
exponent x is, in principle, a characteristic of the universality class the growth process belongs to. 

The divergence of {w2)l suggests that the width (or roughness) distribution PL{'W2)dw2, deflned as the probability 
that W2 is in the interval [w2, W2 + d'W2], is a natural choice when searching for nontrivial distributions. Indeed, if the 
picture about criticality is correct then (102) l gives the only relevant scale in the problem and thus it follows from 
dimensional analysis that Pl{w2) can be expressed in the following form 

Pl{w2)^-^'^(-^] (2) 

{W2)^ \{W2)r.J 
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where ^(x) is a universal scaling function, the object of our main interest in this talk. 

From a utilitarian point of view, three questions should be immediately answered. First, can this quantity be 
measured in experiments; second, can we calculate it theoretically; and third, are the <i>(x)-s of different universality 
classes sufficiently different to be distinguishable. The answer to all three questions is a yes. Present day experiments 
can measure surfaces at high resolution J30J and thus the P{w2) distribution can be built As to the theoretical 
calculation, if a model can be simulated then P{w2) can, of course, be measured. In simple cases, however, one may 
know the nonequilibrium steady-state distribution V[h{r)\ ^ exp{— S'[/i]} and $ can be calculated exactly (an example 
will be presented in Sec lIII C|) . Here we sketch only the first steps of the calculation. 

Formally, P{w2) is obtained from V[h{f)\ as a path integral 



P{W2) = / Vh{r)5{w2 



h ]]ex^]){~S[h]} 



(3) 



where the overbar ft," denotes spatial averaging. In practice, it is easier to calculate the generating function 



G{s) 



P{w2)dw2=JV / Vh{f) exp{-S[h]- s[h^ ~ h ]} , 



(4) 



where A/" is a normalization constant. The above expression is instrumental in finding models where P(w2) can 
be obtained analytically. The path integral in Eq.Q) is the partition function of a model with an effective action 

Seff[h] = S[h] + s[h'^ ~ h ]. The terms added to S[h] are quadratic functionals of h and so one expects that the 
generating function and thus P{w2) can be evaluated exactly if the original model defined by S{h] is solvable. There 
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FIG. 2; Fluctuation distributions for the d = 1 dimensional Edwards- Wilkinson (EW) and MuUins-Herring (MH) models, and 
for the d = 2 EW and Kardar-Parisi-Zhang (KPZ) models. The scaled probability of the roughness $ = (102)^(^2) is plotted 
against the scaled roughness x = 1112 /{w2). The d — 1 scaling functions are the limiting functions when the size of the system 
L goes to infinity. For the d = 2 EW model, a large but finite-L scaling function is plotted since one obtains a delta function 
in the L ^ 00 limit. The KPZ distributions were built by simulating L x L systems. 



are several growth models which are exactly solvable since S[h] is a quadratic functional of h. A notable example 
is the Edwards- Wilkinson (EW) model describing growth governed by surface tension and having S[h] ^ (Vft.)^. 
Another example is the MuUins-Herring (MH) equation 8] modeling curvature driven growth and having a steady 
state characterized by S[h] ^ (A/i)^. These models belong to distinct universality classes since (w2)l ^ L^^ with 
XEW,d=i = 1/2 while XMH,d=i = 1 in d = 1 dimension. Accordingly, their scaling functions should be different. The 
d = \ scaling functions for these models 17, 19, 23;! are displayed on Fig[21and one can see that they are easily 
distinguishable . 

Fig|21 also shows the <I>-s for the d — 2 EW and KPZ models [l^ I2M |26| (the latter is a nonlinear model dis- 
cussed in Sec lIII B)) which are again in different universality classes since XEW,d=2 = (logarithmic divergence) while 
XKPZ,d=2 ~ 0.39. As can be seen from Fig|3the scaling functions strongly differ from each other and, furthermore, 
they are also distinct from the $-s of the d — 1 models. Thus Fig|21 gives an "answer by example" to the third 
utilitarian question posed above. 

We conclude this section with two notes on Fig|21 First, it should be remarked that the d = 2 EW scaling function 
approaches a delta function on the scale of 1^2 /(W2), meaning that the fluctuations of W2 do not diverge in the L ^ 00 
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limit. The delta function, however, hides an interesting structure which can be revealed by an appropriate choice of 
the scaling variables Second, the collapse of the L = 11 and L = 157 KPZ results demonstrates an important 
point about the explicite L-dependence of the <i>-s. Indeed, the scaling functions have finite size corrections, i.e. they 
have explicite i-dependence in addition to the L-dependence through the argument x — W2/ {w2)- Not much is known 
about the finite-size correction of <l>-s but the experience with a large number of models indicate that the explicite 
L-dependence is negligible when the number of surface sites becomes smaller than the number of sites in the bulk. 
This is what we see in Fig|21on the example of the d = 2 KPZ model. 



III. APPLICATIONS 



A. Discovering structures and similarities in experiments and in simulations 

A straightforward way of applying the picture gallery built for surface models is to take the results of a surface 
growth experiment, build the scaling function of the width distribution, and compare it with the pictures in the 
gallery. If one finds agreement with one of the pictures then one can reason about the physical processes which are 
relevant in the given growth process. These types of arguments can be found for example in Ref . . 

A more sophisticated application was the establishment of a connection between the dissipation fluctuations in a 
turbulence experiment and the magnetization fluctuations in the d = 2 XY model at low temperatures jsH . Since 
the low-temperature fluctuations in the d = 2 XY model are equivalent to the surface fluctuations d = 2 Edwards- 
Wilkinson model, the discovery of the above connection prompted a search for an interface interpretation of the 
dissipative structures in the turbulent system j32l |. 

In more theoretical applications, the scaling functions developed for surfaces have been used to find the univer- 
sality class of massively parallel algorithms and thus to establish their scalability ^3§\. Furthermore, they were also 
instrumental in establishing the universality class of fronts propagating into unstable phases [s^ l . Below, we describe 
in detail two more theoretical applications. The first is notable for its logics of approaching a controversial problem 
while the second is remarkable for the puzzle in the end result. 



B. Upper critical dimension of the Kardar-Parisi-Zhang model 

A nontrivial feature of the approach advocated in this talk is that the building of the PDF-s does not involve any 
fitting or other procedures involving subjective judgment. In order to appreciate the advantages of this feature, we 
shall discuss below theproblem of upper critical dimension of the KPZ equation. 

The KPZ equation 35] is the simplest nonlinear generalization of the EW model. In addition to the surface-tension 
effects, it also takes into account that the surface grows along its normal provided the attachment dynamics is isotropic. 
Then, as one can see from Fig^ the vertical velocity of the surface has a correction term proportional to (Vh)'^, and 
the equation, in lowest order in the nonlinearities, can be written as 

dth = vV^h + X(yh)^ +7]. (5) 

Here v and A are parameters, while ri{r, t) is a Gaussian white noise. The above equation has been investigated inten- 
sively since it gives account of a number of interesting phenomena (Burgers turbulence, directed polymers in random 
media, etc.). Nevertheless, a number of issues remained unsolved. In particular, there is no agreement on upper 
critical dimension above which a mean-field theory would be valid. Mode-coupling and other phenomenological 
theories (see Ref.[36l| and references therein) suggest that = 4 while all the numerical work (see e.g. Ref. (33) fail 
to find a finite du- There are, of course, problems with the approximations in the phenomenological theories as well as 
with the multiparameter fits of the simulation results, and the debate is rather controversial. We show now, following 
the lines of our recent work [2^ , how the scaling functions of the roughness can shed some light on this controversy. 

The logics of the solution originates in the theory of critical phenomena. It is known that the universal scaling 
functions depend on d for d < du while their dependence is lost in the mean-field region {d > du)- Thus finding 
smoothly varying <I>-s in dimensions d — l,d,d + 1 should lead to the conclusion that d du- To use this argument, 
we compare on FiglSlthe d = 3, 4, 5 KPZ scaling functions using simulation data for restricted solid-on-solid (RSOS) 
growth models 26, 37]. The message from Figgis clear: It is highly unlikely that du = 4 would be the upper critical 
dimension of the KPZ equation. Note that we arrived at this result without using any approximation or fitting 
procedure. The only way out of the conclusion d ^ du would be if there would be some finite-size corrections to the 
scaling functions which persist at large sizes. One cannot see such corrections for d = 2 in Fig|21and the finite-size 
analysis of the 'I'-s in 3 < d < 5 leads to a similar conclusion (2^ . 
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FIG. 3: Scaling functions for the roughness distributions in the d = 3 — 5 dimensional KPZ models. The scaled probability of 
the roughness $ — (iU2)f (^2) is plotted against x — W2/{w2). The characteristic linear size of the simulated systems is given 
by L. 



C. 1// noise and the Fisher-Tippett-Gumbel distribution 

Studies of the scaling functions often yield connections which are intriguing but hard to understand. The aim 
of this section is to derive an example of such puzzling result. Namely, we shall show that the PDF of the mean- 
square fluctuations in 1// noise is related to the Fisher-Tippett-Gumbel distribution which is one of the limiting 
distributions of extreme value statistics. The derivation in itself is worth going through since it demonstrates the 
steps of calculating P{w2) in a case where the usual choice of the scaling variable ^2/(^2) would yield a delta- function 
PDF in the large-size limit. The derivation is an expanded version of what can be found in a recent letter of ours '20|. 

1// noise is usually associated with a time signal so we should start by imagining the interface in FigQ as a time 
signal i.e. we should make the replacements a; — > i and L T with T being the period of the signal. The mean 

2 

square fluctuations W2{h) = [h{t) — h] are given by averaging over period [0,r]. Then the expressions for P{w2) 
and for the generating function G(s) [Eqs.© and are unchanged and one is left with the problem of finding a 
suitable S[h\ representing the 1// noise. In Ref . . we proposed that the path probability of a Gaussian, periodic, 
random phase, and perfectly 1// noise with the dispersion being linear for all frequencies, can be described by the 
following action 

N N 

S = cy |«||c„P = 2a^n|c„|2, (6) 

n=-N n=l 

where a is an effective surface tension in the language of surfaces and the c„-s are the Fourier coefficients of the 
periodic [ h{t + T) = h{t) ] signal 

N 

h{t)= c„e2""*/^, c_„ = < . (7) 

n=-N 

Several notes are in order here. First, there is a cutoff in above Fourier series. Its meaning is that the time signal is 
resolved at time increments of t = T /N . Using N finite makes the steps of the calculation simple and the N 00 
(equivalent to the T — s- cx) limit) can be taken in the final results. Second, the power spectrum with the above action 
is indeed 1// 

(|c„P)^l/|n|, (8) 
and, third, there is a simple meaning to W2 in the noise terminology since 

N 

W2^2Y\CnW (9) 
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i.e. W2 is the integrated power spectrum. 

Turning now to the evaluation of P{w2), let us note that the functional integral I^J can be written in terms of the 
Fourier amplitudes as 



G(s) - AT J] y dc„ J del exp 



N 

^ 2{am 

m— 1 



(10) 



where A/" is a constant to be determined from the normalization {G{Q) — 1) condition. Carrying out the integrals, one 
finds the generating function in terms of a product 



N 

G{s) = n ^ • (11) 



n=l 



The structure of the above expression is common to various Gaussian growth models where S[}i\ is a quadratic 
functional of h. For example, to obtain the G{s) for the EW model, one should just make the substitution n ^ n? 
and (T — > a/N . Since G{s) has poles on the negative real axis, there is a straighforward method for calculating P(w2)- 
Namely, the inverse Laplace transform of G(s) is an integral along the imaginary axis, 



IOC 



-too "=1 

and thus P{w2) is obtained as a sum of contributions from the poles. The method works p^ll7ll2^ but in our case it 
yields a rather complicated expression which can be shown to converge to a delta function if the usual scaling variable 
X = W2/ {W2) is used. Actually, this result can be understood by just calculating the first two cumulants of W2. The 
first cumulant diverges for — > 00 as 

1 ^ 1 

-^n-i w - [lniV + 7] (13) 



dC 

{W2) ■ 



ds 

where 7 = 0.577... is the Euler constant, while the fluctuations of W2 are finite 



{wi)-{w2?^y^; a=^. (14) 

This means that if we scale W2 with (^2) then the width of the distribution goes to zero for N 00 hence the 
conclusion about the delta function. In order not to lose information about the possible structure at small 11)2 — {VJ2), 
one should introduce a scaling variable which expands the delta function. This can be achieved by introducing the 
following variable 

W2_^^_ (15) 

Substituting W2 from the above expression into Ea. (|12|l and using Eas. (|13|l and (|14|l . one finds that the limit N 00 
can be taken and a scaling function in terms of y emerges 



zoo 



<i>(y) ^ ^{wj) - {W2)^P{W2) = / l^e^n^^. (16) 

— too 

Noting that infinite product above can be expressed [s^ through the T function 

00 _s_ 

n^^.-/"r(i + i), (17, 



71=1 

and using Euler's integral representation for the F function, one finds a surprisingly simple analytical result 

$(y) = ae-(''^+'')~'="'°''^^' . (18) 

The real surprise actually is that $(j/), shown in Fig. ^ is the the so called Fisher-Tippett-Gumbel distribution [s^ 
which is one of the three limiting forms of extreme value statistics. 

It should be emphasized that we do not see any physical reason why should there be a connection between the 
fiuctuations of a 1// signal and the extreme value distributions. It seems to be a puzzle, however, whose solution 
may reveal an interesting picture. 
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FIG. 4: The Fisher- Tippett-Gumbel distribution. The inset shows it on semilogarithmic scale to demonstrate the exponential 
decay at large arguments. 

IV. EFFECT OF BOUNDARY CONDITIONS 

As a final topic, I would like to consider the boundary-condition dependence of the scaling functions. As discussed 
in previous sections, the universality is due to the diverging fluctuations which in turn imply that the correlation 
length is infinite in the system. The infinite correlations mean that the boundaries are felt in the bulk and thus it 
is not entirely unexpected that the scaling functions are sensitive to the boundaries. The important question is how 
large the boundary effects are. 

The effects of boundaries have been studied 0| for Gaussian signals with 1//" power spectrum. It was found 
that the boundary effects are large for a > 4, easily noticeable in the range 1 < a < 4, and they disappear entirely 
for a < 0.5. In order to demonstrate the magnitude of the boundary effect in the physically most relevant range of 
1 < a < 4, I would like to consider the problem of the randomness of the digits of tt which, as we shall see, corresponds 
to an a = 2 problem. 

The statistical properties of the digits of tt appears to occupy the mind of a few mathematicians "ZS] ^ . Although 
no rigorous proof exists yet, it is believed that the digits are random. This belief is based on generating and analyzing 
a large number (10^ — 10*) of digits. We shall test the belief of randomness by mapping the digits onto an interface 
(see FiglS)) and calculating the scaling function of the width of the interface. 



31415926535! 
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FIG. 5: Mapping between the digits of tt and a configuration of the so called rooftop model |4lll42| of interface. The slope of 
the interface can be only ±1 with —1 occurring for digits — 4 and -1-1 for 5 — 9. 

If the digits of tt are random then the interface shown in Figgis just a random walk equivalent to the steady state 
of the Edwards- Wilkinson model as well as to Gaussian signals with 1/ p power spectrum. For periodic 1/ p signals 



^ See e.g. |http://www.sfu.ca/~ pborwein| 
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the width distribution was calculated in Ref. [1^]. It is not straightforward, however, to make a comparison with 
this distribution function since the digits of tt provide only a long but single signal (not unlike to some experimental 
situations). A histogram of roughness can be built, nevertheless, by calculating W2 for segments (windows) of the 
signal whose size is much smaller than the total length but at the same time large enough that the finite-size effects 
in the scaling function would be small. Clearly, the boundary conditions for the segments are not periodic (we call 
them window boundary conditions, WBC). The PDF of the roughness using WBC turns out to be a well defined 
function in the limit of the total length of the signal going to infinity, and, furthermore, this function is found ^2l] to 
be independent of the boundary conditions for the full signal. Comparison of this function with the width distribution 
of TT is displayed on Fig|Hl where one can also see the scaling function for periodic boundary conditions (PBC). The 
difference between the PBC and WBC scaling functions is easily observable and it is clear that we would have 
concluded incorrectly about the non randomness of the digits of tt had we ignored the boundary condition dependence 
of the scaling functions. 




FIG. 6: Distribution of the roughness of tt calculated by mapping its first 10^ digits onto the rooftop interface model (Figl^. 
The scaled distribution $ W2 > P{uj2) is obtained by using a measuring window of 10^ digits and the scaled variable is 
X = 1112/ < W2 >■ The comparison is with random walks with periodic- and window (free) boundary conditions. The collapse 
with the window BC curve is in agreement with the notion that the digits of n are random. 

The moral from the story of tt is twofold. On one hand it cautions us that the scaling functions are useful theoretical 
instruments only if the boundary conditions for experimental data are carefully specified both at the measuring and 
at the analyzing stage. On the other hand, the boundary-condition dependence of the scaling functions suggests that 
they "/tear the shape of the drurri^ i.e. they may be used to see the shape of objects which are embedded in other 
media and can be seen only through their fluctuations. 



V. FINAL REMARKS 



It should be clear that only a small part of the picture gallery was exhibited in this lecture. The gallery is far from 
complete and we expect that scaling functions originating from newly discovered nonequilibrium universality classes 
will enrich it regularly. Furthermore, we expect that limiting distributions such as the ones emerging in extreme 
statistics will also be included and will have much wider use in physics. Finally, strongly fluctuating quantum systems 
may also provide valuable novelties. 
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